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In this work, we calculate the branching ratios, polarization fractions and CP asymmetry pa- 
rameters of decay modes B p{uj)p{uj) in the perturbative QCD approach, which is based on 
kT factorization. After calculation, we find the branching ratios of B'^ — > p~ , B^ P'^ and 

B^ — > p+tj are at the order of 10~^, and their longitudinal polarization fractions are more than 
90%. The above results agree with BaBar's measurements. We also calculate the branching ratios 
\^ • and polarization fractions of p'^p'^, B° p^uj and B^ — > ujui decays. We find that their 

' longitudinal polarization fractions are suppressed to 60-80% due to a small color suppressed tree 

, contribution. The dominant penguin and non-factorization tree contributions equally contribute to 

' the longitudinal and transverse polarization, which will be tested in the future experiments. We 

^ ' predict the CP asymmetry of i?" — > P^ P^ and B^ — > p'^ fP , which will be measured in B factories. 
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I. INTRODUCTION 



> . 

' The study of exclusive non-leptonic weak decays of B mesons provides not only good opportunities for 

testing the Standard Model (SM) but also powerful means for probing different new physics scenarios beyond 
the SM. The mechanism of two body B decay is still not quite clear, although many physicists devote to 
I this field. The hadronic effects must be important while a reliable calculation of these effects is very difficult 
l/^ . Q . Starting from factorization hypothesis Q , many approaches have been built to explain the existing data 
' and made some progress such as generalized factorization Q, QCD factorization (BBNS) H, perturbative 
QCD Approach (PQCD) |^, and soft-coUinear effective theory (SCET) j^. These approaches separately 
' pi explained many of the B PP and B PV decays though some flaws existed in different approaches. 
I I Recently, B — > VV decays such as i? ^ (j)K* dj, i? ^ pK* j^], have aroused many interests of physicists. 

It is known that both longitudinal and transverse polarization states are possible in i? ^ VV decay modes. 
(— I I So, the theoretical analysis of i? ^ VV is more complicated than B PP and B PV. The predictions of 
those decays' polarization fractions according to the naive factorization do not agree with the experimental 
results, although many ideas l^ llol[ll[ have been proposed to explain this phenomenon. Some people think 
, that it is a signal of new physics [l^illdl- Very recently, both BaBar and Belle have measured the branching 
^ • ratios and polarizations of the decays B^ P^P~ and i?+ P^P^ ^Mi^M>^M>^^j some decay modes have 
very large branching ratios. The longitudinal polarization fractions are also very large, which are different 
from that of B ^ (j)K*. 

In this paper, we will study the branching ratios, polarization fractions and CP violation parameters of 
B p{uj)p{uj) decays in the PQCD approach. At the rest frame of B meson, the B meson decays to light 
vector mesons with large momentum. Because the two light mesons move fast back to back, they have small 
chance to exchange soft particles, therefore the soft final state interaction may not be important. A hard gluon 
emitted from the four quark operator kicks the light slow spectator quark in B meson with large momentum 
transfer to form a fast moving final state meson. Therefore, the short distance hard process dominates 
this decay amplitude. In this factorization theorem, decay amplitude is written as the convolution of the 
corresponding hard parts with universal meson distribution amplitudes, which describe non-perturbative 
hadronic process of the decay. Because the Sudakov effect from fc^ and threshold resummation 18], the end 
point singularities do not appear. 

This paper is organized as follows. In Sectional we give some ingredients of the basic formalism. The 
numerical results for branching ratios and CP asymmetry are given in Section IIIII and IIVI respectively. We 
summarize our work at Section [V] 
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II. FORMALISM 



The recently developed PQCD approach is based on kx factorization scheme, where three energy scales 
are involved [3 ■ The hard dynamics is characterized by yJmB^QCD, which is to be perturbatively calculated 
in PQCD. The harder dynamics is from mw scale to scale described by renormalization group equation 
for the four quark operators. The dynamics below ^Jms^QCD is soft, which is described by the meson wave 
functions. The soft dynamics is not perturbative but universal for all channels. Based on this factorization, 
the B ^ pp decay amplitude is written as the following factorizing formula |l9j |. 

A4 ^ Jdxidx2dx:}bidbib2db2b:}db^ 

xTr[C(i)$B(a;i, 5i)$p(x2, 62)$p(x3, 63)^^(2:., 5„ i)^t(x,)e-^(*)] , (1) 

where Tr denotes the trace over Dirac and color indices. C{t) is Wilson coefficient of the four quark 
operator which results from the radiative corrections at short distance. The wave function $m absorbs 
non-perturbative dynamics of the process, which is process independent. The hard part H is rather process- 
dependent and can be calculated in perturbative approach, t is chosen as the largest energy scale in the 
hard part, to kill the largest logarithm. The jet function St{xi), called threshold resummation, comes from 
the resummation of the double logarithms In^x^. The Sudakov form factor S(t), is from the resummation 
of double logarithms In^ Qb 5^, 22} ■ 



A. Wave Function 



In this paper, we use the light-cone coordinates to describe the four dimension momentimr as ,P^)- 
The B meson is treated as a heavy-light system, whose wave function is defined as: 

,(in) 



mb0j{Q)d^,{z)\B{p)) 
iS,,i 



dxd kxe' 



i{xp —hTZj 



' + MB)lz(j)B{x,kT)]^p , 



(2) 



where the indices a and (3 are spin indices, i and j are color indices, and Nc = 3 is the color factor. The 
distribution amplitude (ps is normalized as 



dxi(j)B{xi,bi = 0) 



ft 



(3) 



where bi is the conjugate space coordinate of transverse momentum fcy, and Jb is the decay constant of the 
B meson. In this study, we use the model function 



j(a;, 6) = A^bx^(1 - x)^ exp 



1 f xKlE 
^yujB 



(4) 



where Nb is the normalization constant. We use ujb = 0.4 GeV, which is determined by the calculation of 
form factors and other well known decay modes 

As a light-light system. The p^ meson wave function of the longitudinal part is given by 



P ,al3.ij 



{pip,eL)\dp,{z)uUm) 



(5) 



c Jo 



The first term in the above equation is the leading twist wave function (twist-2), while the others are sub- 
leading twist (twist-3) wave functions. The p meson can also be transversely polarized, and its wave function 
is then 



<p {p,eT)\dpj{z)uai{Q)\0 > = 



dxe"P-^{^T[M'^p{x)+mp^;{x)] 



p ■ n 



(6) 



3 



where n is the moving direction of p particle. Here the leading twist wave function for the transversely 
polarized p meson is the first term which is proportional to 

The distribution amplitudes of p meson, (pp, , 0^, (pj , (p^, and 0°, are calculated using light-cone QCD 
sum rule |2Q||: 



3/p 



■.x{\~x) [l + O.lSCf ^(2x - 1) 



■I p 



{i{2x - if + 0.3(2a; - lf[b{2x - 1)^ - 3] 



+0.21[3-30(2j;-l)2 + 35(2a;-l)'*]} , 
3f^ 

— ^(1 - 2x) [1 + 0.76(10x2 - lOx + 1) 



2J2Nr 



3/, 



" x{\-x) l + Q.2Cl'^{2x~\) 
3 



/2N, 
fp 



[1 + (2x - ly] + 0.24[3(2x - 1)^ - 1] 



2y^2Nr. [ 4 
+0.12[3 - 30(2x - if + 35(2x - 1)^] 



3/p 



AV2Nr 



(1 - 2x) [l + 0.93(10x2 - lOx + 1)] 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



with the Gegenbauer polynomials, 



i(3i2-l), cy2(t) = i(35t4-30i2 + 3), 
^(5^2-1), Cf2(t) = ^(21^4 -14^2 + 1) 



(13) 



B. Perturbative calculations 

For decay B pp, the related effective Hamiltonian is given by |23| 

i/cff = ^ [vudV:, [Ci{p)0^{p) + C2{p)02{p)] - V:,VtdY.C..{p)0,{p)\ , (14) 



i=3 



where Ci{p){i = 1, • • • , 10) are Wilson coefficients at the renormalization scale /x and the four quark operators 
Oi(i = 1, • • • , 10) are 

Oi = (JhUj)v -A{ujdi)v -A, O2 = (hiUi)v~A{ujdj)v~A, 

O3 = {hdi)v-A Y.q%1])v-A, O4 = {bidj)v-A Y.q%<li)v-A, 

O5 = {hd,)v-A T.q%1o)v+A, Oe = {b,dj)v-A T,qi^3li)v+A, (15) 

O7 = ^{hd,)v-Aj2q(^qi^jl3)v+A, Os = f (Mj ) V-A Eg 6,(9,9^) V+yl, 
Og = ^{bidi)v-Aj2q(^qi^jl3)v~A, Oio = ^{bidj)v~Aj2q(^<]i^j1i)v-A- 

Here i and j are SU{3) color indices; the sum over q runs over the quark fields that are active at the scale 
p. = Oirrii,), i.e., q £ {u, d, s, c, 6}. Operators Oi,02 come from tree level interaction, while 03,04,05,05 
are QCD-penguin operators and 07,08,09, Oio come from electroweak-penguins. 

Similar to the B — > tttt decays pj, there are eight types of Feynman diagrams contributing to i? ^ P~ 
decay mode at leading order, which are shown in Fig.^ They involve two types: the emission and annihilation 
topologies. Each type is classified into factorizable diagrams, where hard gluon connects the quarks in the 
same meson, and non-factorizable diagrams, where hard gluon attaches the quarks in two different mesons. 
Through calculating these diagrams, we can get the amplitudes Mj^, where H = L^N^T standing for the 
longitudinal and two transverse polarizations. Because these diagrams are the same as those of i? ^ (pK* 
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and S ^ K*K* the formulas oi B ^ pp are similar to those of B ^X* or B ^ K*K*. We just 
need to replace some corresponding wave functions, Wilson coefficients and corresponding parameters. So 
we don't present the detailed formulas in this paper. 




P 



FIG. 1; The leading order Feynman diagrams for B pp. 



III. NUMERICAL RESULTS FOR BRANCHING RATIOS AND POLARIZATIONS 



In our calculation, sortie parameters such as meson mass, decay constants, the CKM matrix elements and 
the lifetime of B meson are given in Table 

Taking B'^ P^P^ as an example, we know that {H = L, N, T): 



Mh - V:.VudTH-VXVtdPH 



(16) 



with definition: CKM phase angle a — arg 



and zh 



Ph 
Th 



The strong phase 5^ and 



ratio Zh of tree (T) and penguin [P) are calculated in PQCD approach. In PQCD approach, the strong 
phases come from the non-factorizable diagrams and annihilation type diagrams because quarks and gluons 
can be on mass shell. Numerical analysis also shows that the main contribution to the relative strong phase 
5h comes from the penguin annihilation diagrams. B meson annihilates into qq quark pair and then decays 
to pp final states T, "24"| . In hadronic picture, the intermediate qq quark pair represents a number of resonance 
states, which implies final state interaction. These diagrams also make the contribution of penguin diagrams 
more important than previously expected. 
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TABLE I: Parameters we used in numerical calculation [2^ 



Mass 


mgo = 5.28GeV = 5.28GeV 
nip = 0.77GeV = 0.78GeV 


Decay 
Constants 


/s = 196 MeV fp = U = 200MeV 
= /i = leOMeV 


CKM 


IKd! = 0.9745 \V^b\ = 0.042 
= 0.0025 \Vtb\ = 0.999 


Lifetime 


Tgo = 1.54 X 10"^^ s tb+ = 1.67 X lO^^^s 



In Table ITTl we show the numerical results of each diagram in 5° P^P^ decay mode. From this table, 
we find that the most important contribution (about 95%) comes from the first two factorizable emission 
diagrams Fign(a) and (b), especially for the longitudinal part. But the first two diagrams can not contribute 
to the relative strong phases. The main source of strong phases are from the annihilation diagrams, especially 
penguin diagrams of Fign(e)(f) and (g)(h). We can calculate that the strong phase for each polarization is 
6l = 13.6°, Sn = 42°, and St = 39°. 

TABLE IL Polarization amplitudes (10^'^GeV) of different diagrams in B° P^P~ decay 



Decay mode 


(a) and (b) 


(c) and (d) 


(e) and (f) 


(g) and (h) 


L(T) 


77 


-2.4 + 0.6j 





-1.4 - 3.4i 


L{P) 


-3.1 


0.14 + 0.03i 


3.0 - 1.7i 


0.39 + 0.57i 


N(T) 


8.7 


1.3-0.05i 





0.04 - 0.09i 


N{P) 


-0.34 


-0.03 + 0.007i 


1.6 + 0.8i 


-0.002 + 0.009i 


T(T) 


17 


2.7 - 0.004i 


0.04 + O.Olj 


0.002 - 0.008i 


T{P) 


-1.8 


-0.07 + 0.02i 


3.2 + 1.7i 


0.004 + 0.004i 



In the same way, we can get the formula for the charged conjugate decay B'^ p : 

Mh = VubV:aTH{l + ZHe'("+*«)). (17) 

Therefore, the averaged branching ratio for B p^ p^ is: 

Ml ^ \V:f,VudTH\^ {1 + 2zH cos a cos Sn + zj,). (18) 

Here, we notice the branching ratio is a function of cos a. This cos a behavior of the branching ratio is shown 
in Fig[21 In principle, we can determine angle a through ea. H18|) . However, the uncertainty of theory is so 
large (also shown in Figl^J to make it unrealistic. Firstly, the major uncertainty comes from higher order 
correction. In the calculation of S ^ Kit |2^, the results show that the next-to-leading order contribution 
can give about 15% — 20% correction to leading order. Secondly, the wave functions which describe the 
hadronic process of the meson are not known precisely, especially for the heavy B meson. Using the existing 
data of other channels such as i? — > irli/ 26], i? — > Dtt 27], B — > Ktt, tttt 5], etc, we can fit the B meson wave 
function parameter lub = 0.4 ± 0.1. Another uncertainty comes from parameter c of threshold resummation 
[34], and it varies from 0.3 to 0.4. In leading order, considering the uncertainty taken by lu and c, we give the 
branching ratios and polarization fractions in Table IhH together with averaging experimental measurements 

HI El 111 Ea. 

There are still many other parameters existed such as decay constants, CKM elements, and we will not 
discuss the uncertainty here. The polarization fractions of these decay modes are not sensitive to the above 
parameters, because they mainly give an overall change of all polarization amplitudes, not to the individ- 
ual noes. From our calculation, we find that these polarization fractions are sensitive to the distribution 
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amplitudes of vector mesons. However, the distribution amplitudes we used are results from Light-Cone 
Sum Rules [20l|. which are difficult to change. Anyway, 20% uncertainty from the meson distribution ampli- 
tudes for the polarization fractions are possible. The range of CKM angle a has been well constrained as 
a = (98^5 5)° 2§j, so that its small uncertainty affects very little on the branching ratios. 




25 50 ■75 Ijj 125 150 175 

a(degree) 



FIG. 2: Average branching ratio with theoretical uncertainty of B p as a function of CKM angle a, where 

the shaded band shows the 1 cr constraint for a 



TABLE III: Branching ratios and polarizations fractions of B — > p{uj)p{uj) decays from theory and experiments 
J^|17|1. In our results, the uncertainties come from ujb and c respectively. 





BR (10" 










Channel 


Theory 


Exp. 


Theory 


Exp. 


/||(%) 






35 ± 5 ± 4 


30 ±6 


94 


96j:* 


3 


3 




17±2± 1 




94 


99 ±5 


4 


2 


B+ ~+ p+LO 


19±2± 1 




97 


88t\l 


1.5 


1.5 


B" ^ pOpO 


0.9 ±0.1 ±0.1 


< 1.1 


60 




22 


18 


B» ^ P°LJ 


1.9 ±0.2 ±0.2 


< 3.3 


87 




6.5 


6.5 


B° ^ULJ 


1.2 ±0.2 ±0.2 


< 19 


82 




9 


9 



From above results and Table UTTl some discussions are in order: 

• For simplicity, we set that the have same mass, decay constant and distribution amplitude. In 
quark model, the p° meson is {uu — dd) , while w is {uu + dd) . The difference comes from the 

sign of dd, which only appears in penguin operators, so their difference should be relatively small. 

• For the tree dominant decays, most of the contribution to branching ratio comes from factorizable 
spectator diagram (a) and (b), which are the diagrams contribute to the i? — s- p form factor. For 
example in decay mode P i the dominant Wilson coefficients are C2 ± Ci/3 (order of 1) in 
tree level, which is supported by numerical results. The decay — s- p^uj and _B+ have the 
similar situation. Their branching ratios are all at the order 10~^. 

• For decay B p^ p^ , the Wilson coefficient is Ci ± C2/3 in tree level, which is color suppressed. 
In this work, we only calculate the leading order diagrams, and did not calculate the higher order 
corrections. So, the Wilson coefficients we used are leading order results in order to keep consistency. 
In leading order, the sign of C2 is positive while the sign of Ci is negative, which can cancel each other 
mostly. Thus the branching ratio of _B ^ is rather small. If considering next to leading order 
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corrections, the sign of Ci + C2/3 may change to positive, so the branching ratio may become larger. 
This decay should be more sensitive to next leading order contribution. This is similar to the argument 
of B'^ 7r°7r° decay and 5° p'^uj, ujlu. 

• Comparing our results with experiments (world average), we find both branching ratios and polariza- 
tions agree well with only one exception: B~^ p~^ p^. In fact, this is due to a large branching fraction 
measured by Belle p^. 

BR(B+ ^ p+p°) = (31.7 ± 7.ltlj) X 10-^ (19) 

which does not overlap with BaBar's data 0, 0, 0| 

BR(B+ ^ p+pO) = {22.5tli ± 5.8) x IQ-^. (20) 

We are waiting for the consistent results from two experimental groups. As for the color suppressed 
^0 _j. pOpO^ _^ j^o _j. there are only upper limits now, and our results are still below 

the upper limits. 

• In Ref. l30l |. these decay modes have been calculated in QCD factorization approach. For 
B^ p~^ p~ , the branching ratio they predicted is a bit larger than the experimental data, because 
the form factor they used is V^^^ — 0.338. In PQCD approach 31J, this form factor is about 0.318, 
so our results is smaller than theirs. Similar to above decay, our results in decay _B+ — s- p^p" is also 
smaller than results in QCD factorization approach for the same reason. For decay modes i?*^ — > p^p", 
B^ p^Lo and — > louj, our results are much larger than theirs because the annihilation diagrams 
play very important role, and these parts cannot be calculated directly in QCD factorization approach. 

• Both dominant by color enhanced tree contribution, we can see that the branching ratio of B'^ —^ p^ p^ 
is about two times of that of _B+ P^ P^ ■ But in experimental side, the world average results of these 
decays have not so much difference. Neither QCD factorization approach nor naive factorization can 
explain this small difference. The similar situation also appears in the decays B ^ mr Many 
people have tried to explain this puzzle P,!!^- But for the B ^ pp decays, it is stiU early, since the 
very small branching ratio of B^ — > p^ p^ by experiments contradicts with isospin symmetry. We have 
to wait for the experiments. 

• From the Table HTll we know that longitudinal polarization is dominant in decay B^ P^P~i ~^ 
p'^ p^ and B~^ — > p'^uj, which occupies more than 90% contribution, and is consistent with experimental 
data. These results are also consistent with the prediction in naive factorization because the 
transverse parts are suppressed, where Vp — mp/niB- But for B^ — > p^p" decay, the tree emission 
diagrams are mostly cancelled in the Wilson coefficients. As we will see later in Table Hvl the most 
important contributions for this decay are from the non-factorizable tree diagrams in Fig^c) and 
(d) and also the penguin diagrams. With an additional gluon, the transverse polarization in the non- 
factorizable diagrams does not encounter helicity flip suppression. The transverse polarization is at the 
same order as longitudinal polarization, which can also be seen in the column (c) and (d) of Table Hn] 
This scenario is new from the mechanism of the recently penguin dominant process B — > (l)K* |33| . 
where the penguin annihilation guides the dominant transverse contribution. In fact, the iS" —> uju) 
decay has a little larger longitudinal fraction is just due to the fact that there is no non-factorizable 
emission tree contribution for this decay in isospin symmetry. 

Now we turn to discuss the contribution of different diagrams, where — s- p^ p~ and B^ —> p^p^ are taken 
as an example. In the Table Hvl we consider full contribution in line (1), ignore annihilation contribution 
in line (2), without all penguin operator in line (3), and without non- factorization diagram in line (4). 
From this table, we can see that neither annihilation diagrams nor non-factorizable diagrams can change 
the polarization fraction in decay B'^ —>■ p^ p^ ■ They only take about 4% contribution in this decay mode 
just because the emission diagram occupy very large part of the contribution, which can also be seen from 
Table Inl However, the penguin operators especially in annihilation diagrams play an important role in decay 
B" ^ pOp". 

Of course , th e final state interaction is very important in non-leptonic B decays. They can give 0(10"^) 
corrections 0|, but this can not change the branching ratios much for decay modes P^P~ and 

B^ p^ p^ at order 10^^. Thus, in these two decay modes, the final state interaction may not be important. 
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TABLE IV: Contribution from different parts in — > p~ and — > p^p": full contribution in line (1), ignore 
annihilation contribution in line (2), without penguin operators in line (3), and without non- factorization diagrams 
in line (4). 



B° - P+p- 


BR(10-'') 




/||(%) 


/x(%) 


(1) 


35 


94 


3 


3 


(2) 


35 


94 


3 


3 


(3) 


32 


94 


3 


3 


(4) 


38 


96 


2 


2 


B« ^ pOp" 


BR(10-'') 


fL{%) 


/||(%) 


/x(%) 


(1) 


0.94 


60 


22 


18 


(2) 


0.38 


42 


26 


32 


(3) 


0.25 


18 


41 


41 


(4) 


1.18 


83 


8.5 


8.5 



However, in decay i?" — > p''\ the final state interaction may afford larger contribution than our calculation 
(IQ-'^-lO-^), that's to say, our perturbative part may not be the dominant contribution. Although probably 
important, the hadronic effects are not intensively discussed in this paper, since they are beyond the topics 
of our PQCD approach. The contributions of these two sides can be determined by experiments. 



IV. CP VIOLATION IN B° ^ p+ p- AND B+ p+p"(c^) 

Studying CP violation is an important task in B physics. In this section, we discuss the CP violation 
in — > p~ and B^ — > p'^p^{uj) decays. The uncertainty in B^ —> p^ p^ (lu) , ujuj for branching ratios is 
so large that we will not discuss their CP violation here, though it is also very important. In decay modes 
B^ — > p~^ p~ and B~^ —> p'^p^{uj), longitudinal part occupy nearly 95% contribution. So we will neglect the 
transverse parts in the following discussions. 

Using Eqs. (|16I17(I . the direct CP violating parameter is easily derived as a function of CKM angle a. 

\M+\^-\M-f 



^CP 



iM+r + |Af-r 

— 2zsinQ!sin(5i 
1 + 2zL cos a cos 5l - 



(21) 



which is shown in FigEJ The direct CP asymmetry is about (— 10±4)% in decay P^ P - However, the 

direct CP in decay P^ P^ is almost zero, because there is no QCD penguin contribution while elec- 

troweak penguin contribution is rather negligible. On the other hand, because of large penguin contribution, 
the direct CP in B^ p+w is about (—23 ± 7)%, which is even larger than B'^ P^P - The uncertainty 
in above results come from 90° < a < 110° and 0.3 < lub < 0.4 in B meson wave function. 

For the neutral B'^ decays, there is more complication from the — B'^ mixing. The time dependence of 
CP asymmetry is: 

Acp ^ AcpCOs(Amt) +sin(ATOi)a£+e/, (22) 

where Am is the mass difference between the two mass eigenstates of neutral B mesons. The A'}.% is already 
defined in Ea. H21|) . while the mixing-related CP violation parameter is defined as 

_ -2Im(Acp) 
1 + \acp\ 

where 



V,lVt, < f\H,ff\B\ > 
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2 5 50 



.25 150 175 



FIG. 3: Direct CP violation parameter Aqp as a function of a with ujb ~ 0.4. The solid line is for , 



a (degree) 

iction of a wit' 

dot-dashed for and dashed line for p^ui. The shadow part is a band with 90° < a < 110° 

Using Eqs. H16I17|I . we derive as: 



•lia 



1 + ZLe'^^^-°'^ 

1 + ZLe^(^^+°') ' 



(25) 



Thus, the parameter a^+e' is a function of a, if the penguin pollution is very small, a^+e' is about — sin 2a. 
From the function relation of Fig.^ we can see that a^+^' is not exactly equal to — sin2Q;, because of the 
penguin pollution. 



0.5 



03 




-0.5 



; 25 50 125 iso 175 

a (degree) 

FIG. 4: Mixing induced CP violation parameters a<,+j/ of B'^ P^ P~ as a function of CKM angle a with uj = 0.4. 
The shadow part is a band with 90° < a < 110°. 



If we integrate the time variable t of Eq. H22|) , we will get the total CP asymmetry as 

1 



A. 



CP 



1 + x 



A dir 
2^CP 



1 + X- 



(26) 



with X = Am/T ~ 0.771 for the — mixing in SM Through calculating, we notice that the Acp 
is (-10 ± 4)% with uncertainty from 90° < a < 110° and 0.3 < < 0.4. 
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V. SUMMARY 

In this work, we calculate the branching ratios, polarizations and CP asymmetry of B — > p{uj)p{uj) 
decays in perturbative QCD approach based on kx factorization. After calculating all diagrams, including 
non-factorizable diagrams and annihilation diagrams, we found the branching ratios of B'^ ^md 
5+ P^ P^ are at order of O(10~^), and the longitudinal contributions are more than 95%. These results 
agree with the BaBar's data well. Moreover, we also predict the direct CP violation in P^ P~ and 

P^ 1 and mixing CP violation in B'^ P^ P ^ which may be important in extraction for the angle 
a. The longitudinal polarization for B^ p'^p'^, p'^^, are suppressed to 60% — 80% due to the large 
non-factorizable tree contribution to these decays. These results can be tested in B factories in future. 
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